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We investigate theoretically the low-energy physics of semiconductor Majorana wires in the vicinity of a 
magnetic field-driven topological quantum phase transition (TQPT). The local density of states at the end of the 
wire, which is directly related to the differential conductance in the limit of point-contact tunneling, is calculated 
numerically We find that the dependence of the end-of-wire local density of states on the magnetic field is 
nonuniversal and that the signatures associated with the closing of the superconducting gap at the Majorana 
TQPT are essentially invisible within a significant range of experimentally relevant parameters. Our results 
provide a possible explanation for the recent observation of the apparent nonclosure of the gap at the Majorana 
TQPT in semiconductor nanowires. 



PACS numbers: 03.67.Lx, 03.65.Vf, 71.10.Pm 

The recently reported observation (TJ of a zero bias peak 
(ZBP) in conductance measurements on semiconductor (SM) 
nanowires coupled to superconductors (SCs) may represent 
the first experimental evidence of Majorana fermions (MFs), 
which are theoretically predicted to exist in topological super- 
conductors 1 2]. Topological SC states capable of supporting 
MFs can be realized in SM wires with proximity-induced su- 
perconductivity by driving the system through a topological 
quantum phase transition (TQPT) using a suitably directed 
magnetic field (2H3- At the TQPT, the SC gap necessarily 
vanishes [2-6]. The absence of any signature associated with 
the gap closure casts serious doubt on the Majorana fermion 
interpretation of the ZBP in the recent experiment 0. 

In this Letter, we offer a possible explanation for the ob- 
served non-closure of the gap. By solving numerically an ef- 
fective tight-binding model for multiband nanowires with re- 
alistic parameters we show that, in the vicinity of the TQPT, 
the amplitude of the low-energy states near the ends of the 
wire may be orders of magnitude smaller than the amplitudes 
of the localized MFs. Consequently, the contributions of these 
states to the end-of-wire tunneling conductance and LDOS are 
essentially invisible, which results in an apparent non-closure 
of the gap in these quantities. By contrast, the closing of the 
gap mandated by the TQPT is clearly revealed by other quan- 
tities, such as the total density of states (DOS) and the LDOS 
near the middle of the wire. We also show that this non- 
closure of the gap is non-universal, being dependent on the 
behavior of certain low-energy wave functions, and we iden- 
tify the parameter regimes in which signatures associated with 
the closing of the gap are present in the end-of-wire LDOS. 

In Ref. [2 ] Sau et al. proposed that a spin-orbit (SO) cou- 
pled semiconductor (SM) thin film with Zeeman spin splitting 
and proximity induced s-wave superconductivity could be 
used to realize Majorana fermions. For small Zeeman splitting 
T, the semiconductor is in a conventional (proximity-induced) 
superconducting state with no MFs, while for T larger than a 



critical value T c (corresponding to the TQPT where the SC 
gap vanishes), localized MFs [ 6 - 8 1 exist at defects of the 
SC order parameter. In subsequent works GHSl l9lfT4l it was 
shown that in the ID version of this system - the so-called 
'semiconductor Majorana wire', a direct physical realization 
of the Kitaev model 1 8 ] - zero-energy MF states are trapped 
at the wire ends and protected from regular fermionic excita- 
tions by a large mini-gap [ 3 ] ~ E qp ^ IK, where E qp is the 
proximity-induced bulk SC quasiparticle gap. The semicon- 
ductor Majorana wire, which has recently received consider- 
able experimental attention (HEHIIL allows the detection of 
the zero-energy MF as a sharp zero bias peak in local charge 
tunneling measurements (3j |71 QjD at experimentally realistic 
temperatures T < E qp [3]. Here, we study the signatures of 
nonzero low-energy states in local end-of-wire measurements. 

We consider a rectangular SM nanowire with dimensions 
L x ^> L y ~ L z proximity coupled to an s-wave supercon- 
ductor. A realistic model of the nanowire that includes the 
effects induced by proximity to the SC is solved numerically 
for a set of effective parameters corresponding to InSb fol- 
lowing the procedure described in Ref. l20l (see also the sup- 
plemental material for technical details I2T1 ). The wire has 
a cross section L y x L z « 45 nm x 50 nm and is charac- 
terized by a Rashba spin-orbit coupling a = 0.2 eVA. An 
external Zeeman field T = g* /i#£?/2, where B is the mag- 
netic field and g* = 50, is applied along the ^-direction. We 
couple the wire to an s- wave SC with a bulk gap Ao = 1.5 
meV. The proximity-induced effective pair potential in the SM 
is A = 0.25 meV. With increasing T, the wire evolves from 
a non-topological SC state with no MF to a topological SC 
state with MFs localized near the ends via a TQPT at T — T c . 
The SC quasiparticle gap induced in the wire must vanish at 
the TQPT (3l[6l|20|. Such closing of the bulk gap is clearly 
visible in the total DOS, as shown in the top panel of Fig. [T] 

The main finding of this Letter is that, despite the vanish- 
ing of the bulk SC gap at the TQPT, as mandated inflexibly 
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FIG. 1. (Color online) 7b/?: Total DOS as a function of B for \i — 
A/ 2. The closing of the bulk gap is clearly visible at the TQPT, 
B = B c « 0.2 T. Middle: LDOS at the end of the wire as a function 
of B. The strong features associated with A ~ 250 /xeV are only 
weakly dependent on For B > B c a peak associated with the 
Majorana bound state is present at zero energy. The LDOS shows 
no visible signature of the bulk gap closing at the TQPT. Bottom: 
LDOS at the middle of the wire. Note the closure of the gap at the 
TQPT and the absence of the zero-energy Majorana peak. Traces for 
selected values of B are provided in the Supplemental Material (21]. 

by the theory [ 3, 6, 20], this gap closure may in fact not be 
visible in charge conductance experiments that aim to probe 
the end- state MFs. We show this by explicitly calculating 
the end-of-wire LDOS, which is qualitatively related to the 
charge current passed through the end of the nanowire cou- 
pled to a normal lead [ 1 ]. This relation becomes a close corre- 
spondence in the limit of point contact tunneling l22ll . Since 
the LDOS can be related more directly to the physics at the 
microscopic level, we focus on this quantity to shed light on 
the recent experimental results. We consider a ID nanowire 
system with four occupied bands, i.e., four pairs of spin sub- 
bands, (see Fig. 2, top panel) and a chemical potential close 
to the minimum of the top band. This is similar to the exper- 
imental situation in Ref . 1 , where only a few sub-bands are 
thought to be occupied (results for the single-band case are 
provided in the Supplemental Material I2T1D . Specifically, we 
have fi = A/2, where the chemical potential is measured rel- 
ative to the energy of the top occupied band at k x = and 
B = 0. The results are shown in middle panel of Fig. [T] 
A zero-energy peak is clearly visible above the critical mag- 
netic field B c , as well as several strong features with energies 
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FIG. 2. (Color online) Top: Spectrum of an infinite non- 
superconducting wire. The lowest three bands (six sub-bands) cross 
the chemical potential at Fermi wave vectors k 7 ^ > 0. lnm -1 , while 
the top occupied band has k 4 ^ 1 <C O.lnm -1 . The fifth band in unoc- 
cupied. Bottom: BdG spectrum for three different values of B. The 
induced SC pair potential is A = 0.25 meV. The three lower-energy 
bands depend weakly on B and have nearly overlapping contribu- 
tions characterized by six minima at energies ^ 250 fieV. By con- 
trast, the top occupied band is strongly B-dependent (three distinct 
curves at k x < 0.02nm _1 ) and is characterized by a gap that closes 
at the critical field B c ~ 0.2 T. Inset: Relevant portion of the phase 
diagram showing the topological SC phase (yellow/light gray) and 
the topologically trivial phase (white). Fig. [T] corresponds to a cut 
with fi = A/2 (path I), Fig. His for fi = 2 A (path III), while Fig. [5] 
corresponds to B(/jl) — B c (/jl)/2 (path II). 

> A = 250/ieV that depend weakly on B. Note that there is 
no visible signature associated with the closing of the gap at 
the TQPT. This type of behavior is very similar to the experi- 
mentally observed dependence of the end-of-wire differential 
conductance on the magnetic field [Q . By contrast, the LDOS 
calculated at the middle of the wire (Fig. [T] bottom) clearly 
shows the bulk gap closing but no Majorana fermion ZBP for 
B > B c . 

To understand the somewhat unexpected behavior of the 
LDOS in Fig. [TJ we consider a nanowire with a chemical 
potential close to the bottom of the fourth band. This situa- 
tion requires a low magnetic field to drive the system into the 
topological SC state, which is consistent with the experimen- 
tal conditions [ 1 ]. The normal state spectrum of the SM wire 
(for B = 0) is shown in Fig. [2] (top panel). The Bogoliubov- 
de Gennes (BdG) spectra associated with the four occupied 
bands are shown in Fig. [2] (bottom panel) for different values 
of the magnetic field. From the weak dependence on B of the 
BdG minima associated with the three lower energy bands, it 
is clear that the prominent and weakly Independent finite en- 
ergy features observed in the end-of-wire LDOS (Fig. [T] mid- 
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die panel) are associated with these bands. By contrast, the 
contribution to the BdG spectra associated with the top band 
(or the 'Majorana band') depends strongly on B. Note that the 
gap in the top band has a minimum at k x = when B < B c , 
vanishes at B c w 0.2 meV, then reopens for B > B c . Thus, 
the low-energy physics in the vicinity of the TQPT and the gap 
closure at B c is controlled by the Majorana band, while the 
bulk gaps at high Fermi momenta associated with the lower 
bands do not close at the TQPT and depend weakly on B. 

Next, we consider a finite wire with L x — 4.5/im and fo- 
cus on the contributions to the LDOS coming from the BdG 
eigenstates associated with the lower bands. In general, the 
contribution to the LDOS at the end of the wire coming from 
a given state n depends on how fast the amplitude of the cor- 
responding wave function ^ n {x) increases as a function of x 
(note that all wave functions must vanish at the wire ends). For 
wires with confinement energy much larger than A, states near 
the Fermi-level that are associated with the lower-bands have 
a large kinetic energy ^ ^> A, V (see Fig. [5] upper panel), 
which is the difference between the chemical potential fi and 
the bottom of the band n. These lower-band states have large 
Fermi wave- vectors and, consequently, are characterized by 
amplitudes rapidly increasing away from the ends of the wire. 
We identify these states as responsible for the weakly mag- 
netic field dependent feature at w 250 /ieV in Fig. [T] Note 
that the lowest energies of the lower-band states have values 
close to the edge of the bulk gap, which disperses weakly with 
B and does not close at the TQPT (Fig. [2] bottom panel). 

By contrast, states associated with the Majorana band have 
long characteristic wave-lengths and, consequently, their con- 
tributions to the end-of-wire LDOS are strongly suppressed 
whenever these states are delocalized bulk states characterized 
by an envelope with a vanishing amplitude near the ends. In 
general, the lowest energy states associated with the Majorana 
band also contain a localized component characterized by an 
envelope that decays exponentially away from the ends of the 
wire. In the non-topological SC phase, we find that the lowest 
energy states are delocalized (i.e., have a negligible localized 
component) when the chemical potential is below a certain 
crossover value /i c (r) ~ O(A) (blue dots in the inset of Fig. 
[2]) and localized when /i > /i c . This delocalized-localized 
crossover is relatively sharp, being characterized by an energy 
scale of order A. Also, we emphasize that there are multi- 
ple localized-delocalized crossovers characterized by values 
of the chemical potential /i c?n slightly above the minimum of 
each band. The striking difference between the lowest energy 
states on the two sides of the crossover (for n = 4) is illus- 
trated in Fig. [3] (red/dark gray lines). Note that in the topolog- 
ical SC phase (B > B c ) the lowest energy state is always a 
Majorana bound state (yellow/light gray lines in Fig. [3]). 

The localized/delocalized character of the low-energy states 
is directly reflected in the end-of-wire LDOS. If the system is 
in the delocalized regime, /i < /i c , the signature of the low- 
energy states associated with the Majorana band is strongly 
suppressed. In particular, there is no visible signature of the 
gap closing at the TQPT. In this case, the dominant features 
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FIG. 3. (Color online) The lowest energy BdG states associated with 
the Majorana band near the TQPT for different values of \i. The 
red (dark gray) lines correspond to B = 0.9B C , while the yellow 
(light gray) curves are for B — 1.1B C . The insets show the BdG 
spectra for an infinite wire near k x — 0. In the topological SC 
phase (B > B c ), the lowest-energy state is the MF state. In the 
non-topological SC phase (B < B c ), there is a crossover from the 
region with \i > \i c ~ A characterized by localized lowest energy 
states (top panel) to the region with \i < \i c « A where the states 
have mostly extended character (middle and bottom panels). The 
corresponding BdG spectra change from a structure with two min- 
ima (/i > /i c ) to a single minimum at k x = (/x < ji c ). 

originate from states in the lower bands, with energies > A 
that depend weakly on the magnetic field. Note that for /i = 0, 
which corresponds to the minimum of the critical Zeeman 
field T c (/i) = A (see the inset of Fig. [2} and is probably 
very similar to the experimental conditions [1 j, the system is 
in the delocalized region, hence no gap closure at the TQPT 
is visible in the end-of-wire LDOS. However, this behavior 
is not universal. Increasing the chemical potential above fi c 
will move system into the a region characterized by localized 
lowest-energy states. Since the energy of the localized states 
is at the bulk gap edge, the LDOS will contain a feature as- 
sociated with the gap edge that will depend strongly on the 
magnetic field and will reveal the gap closure at the TQPT. 
This situation is illustrated in Fig. [4] We therefore predict that 
measurements probing other regimes of nanowire parameters 
may very well observe gap closing signatures, along with the 
emergence of MFs, depending on the details of the system. 
The qualitatively different behavior of the lowest energy 
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FIG. 4. (Color online) Dependence of the LDOS at the end of the 
wire on B for fi — 2 A. The gap-edge feature stemming from the 
localized end state clearly reveals the gap closing at the TQPT. For 
B > B c the finite energy states are de-localized and, consequently, 
there is no feature showing the re-opening of the gap. Traces for 
selected values of B are provided in the Supplemental Material ['21 1. 

BdG wave functions for /i < \i c ~ A as compared to /i > \i c 
can be correlated with the main features of the corresponding 
BdG spectra. As shown in the inset of Fig. [3] (top panel) the 
BdG spectrum for /i > /i c has two well defined minima at two 
non-zero values of k x . This behavior is similar to that corre- 
sponding to the lower three bands (see Fig. [2| lower panel). 
By contrast, as shown in the insets of the middle and bottom 
panels in Fig. [3] the BdG spectra for /i < /i c have only a sin- 
gle minimum at k x = 0. In the cross-over region, /i ~ /i c , 
the lowest energy BdG state for the finite wire is a linear su- 
perposition of localized and extended states. Because of this 
clear pattern, we find that for \i > \i c the end-of-wire LDOS 
generically reveals the dispersion of the bulk gap in the Ma- 
jorana band via the localized end states at the gap edge, while 
for ii < n c the end-of-wire LDOS does not disperse with B. 
This is further illustrated by the dependence of the LDOS on 
the chemical potential in Fig. [5] Note that the LDOS (bot- 
tom panel) shows features corresponding to the bulk gap edge 
only in the region /i > /i c , while it shows no visible signa- 
tures associated with the gap edge for /i < /i c . We emphasize 
that this behavior is dictated by the localization properties of 
the low-energy wave functions in various regions of the phase 
diagram. These properties are not expected to change quali- 
tatively in the presence of a smooth confining potential l23ll 
finite temperature (24), or weak disorder lf25l (see the Supple- 
mental Material for details ??). 

We establish that the LDOS in the semiconductor Majorana 
wire is a non-universal quantity sensitively dependent on the 
details of the system parameters and that it is entirely possi- 
ble for the LDOS not to manifest any signature of the SC gap 
closing through the TQPT. When the chemical potential is in 
the vicinity of the minimum of an arbitrary band n, we iden- 
tify a sharp localized-delocalized crossover characterized by a 
value of the chemical potential (relative to the bottom of band 
n) ji c ~ A, where A is the proximity induced pair potential. 
For B < B c and \i > /i c , the lowest energy BdG states asso- 
ciated with the top Majorana band are localized near the wire 
ends and contribute significantly to the end-of-wire LDOS. 
Since the energy of these localized states is close to the gap 
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FIG. 5. (Color online) Total DOS (top) and LDOS at the end of the 
wire (bottom) as functions of ji. The dependence on the chemical 
potential is along path II (see inset of Fig. [2]). The gap is smaller 
than the pair potential A = 250 /xeV, as shown by the DOS (top 
panel). In the LDOS, a feature associated with the gap edge becomes 
visible only in the localized region, ji > fi c w 250 /xeV (blue dots 
in the inset of Fig. [5]). 

edge, the LDOS reveals the dispersion of the SC gap with B 
and the eventual closing at the TQPT. By contrast, for /i < ji c 
the lowest energy BdG states associated with the Majorana 
band decay near the wire ends. Consequently, their contribu- 
tion to the end-of-wire LDOS is negligible and the LDOS does 
not reveal the bulk gap closing at the TQPT. This behavior is 
robust against disorder, if it is not too strong to destroy the 
topological phase. In the presence of a disorder potential, the 
low-energy extended states will become localized inside ran- 
dom segments of the wire, but, typically, will still decay near 
the ends (see the supplementary section for details). Also, fi- 
nite temperature will generate broadening effects, but will not 
shift the spectral weight. Since the LDOS at the end of the 
wire is related to dl/dV, which is the experimentally mea- 
sured quantity, we expect that for ji < fi c the bulk gap closing 
at the TQPT should not be seen in tunneling conductance mea- 
surements, although the signature of the MFs would clearly 
show up as a zero-bias-conductance peak. This is consistent 
with recent experiments 1 1 ], which observe no feature associ- 
ated with the bulk gap closing at B = B c , yet a ZBP is visible 
fori? > B c . By contrast, charge conductance experiments 
in the regime /i > /i c , or experiments that directly probe the 
bulk gap 1261 , should reveal the dispersion of the bulk SC gap 
and its eventual closing at B — B c . Experiments probing the 
LDOS at the middle of the wire would reveal signatures of the 
bulk gap closing through the TQPT, but no Majorana ZBP. 

This work is supported by DARPA-MTO, NSF, DARPA- 
QuEST, JQI-NSF-PFC, Harvard Quantum Optics Center, and 
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Note Added: After the original submission of this 
manuscript, a number of theoretical papers I27ti35l ap- 
peared addressing various aspects of the experimental work 
on the possible observation of the Majorana zero mode in 
semiconductor-superconductor hybrid structures. 
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Supplemental Material 

Modeling. We consider a semiconductor nano wire with 
rectangular cross section and dimensions L x ^> L y ~ L z . 
The low-energy physics of the nanowire is described by the 
Hamiltonian 

HsM = H + HsOl = ^ tij^ia ^ ~ C iv Ci(T 

+ y [ c Ud x °y c i - c \+8 y °x c i + h - c -] > (!) 

i,d 

where Hq, which includes the first two terms, describes near- 
est neighbor hopping on a simple cubic lattice with lattice con- 
stant a with ta+s = —to, where S are the nearest-neighbor 
position vectors. In Eq. ([T]) the last term represents the Rashba 
spin-orbit interaction (SOI), c\ is a spinor c\ = (c|^, c|^) with 

c\ a being the electron creation operators with spin a, fi is the 
chemical potential, a is the Rashba coupling constant, and 
& = (cr x ,cry,cr z ) are Pauli matrices. Since the number of 
degrees degrees of freedom in a finite wire is large (of the or- 
der 10 T ), yet Majorana physics is basically controlled by a re- 
duced number of low-energy degrees of freedom (of the order 
10 3 -10 4 ), we project the problem into the low-energy sub- 
space spanned by a certain number of low-energy eigenstates 
of H . Explicitly, the eigenstates of H are 

= ft V^TT sin WTi x ^ (2) 

where n = (n Xl n yi n z ) with 1 < n\ < N\, and Xa is an 
eigenstate of the a z spin operator. The corresponding eigen- 
values are 

/ 7rn x irUy 7rn z \ 

e n = -2t ^cos — +cos — +cos — - 3 j - w , 

(3) 

where the chemical potential fio is calculated from the bot- 
tom of the first band. The low-energy subspace is defined by 
the eigenstates satisfying the condition e n < e max , where the 
cutoff energy e max is of the order 75-100meV. Using this low- 
energy basis, the matrix elements of the SOI Hamiltonian can 



FIG. 6. (Color online) Top: Total DOS as a function of the magnetic 
field for a system with single band occupancy. For B — 0, the chem- 
ical potential is set at the bottom of the first band, i.e., \i — and 
^to P = 1. The closing of the bulk gap is clearly visible at the TQPT, 
B — B c « 0.2 T. Bottom: LDOS at the end of a single channel 
nanowire as a function of B. For B > B c a peak associated with the 
Majorana bound state is present at zero energy. The LDOS shows no 
visible signature of the bulk gap closing at the TQPT. Note the ab- 
sence of any strong feature associated with the induced gap A ~ 250 
/ieV. 

be written explicitly as 



(i>na\Hsol\ll>n'a') <*$n a n' 



(-1) 



n x +n' 



yjerer' 



sin 



N x + 1 



sin 



[x y] 



(4) 



cos 



N x + 1 



COS 



N x + 1 



where the second term in the parentheses is obtained from the 
first term by exchanging the x and y indices. Note that the SOI 
Hamiltonian has the structure i^soi = ^soi + ^soi' wnere 
the first term represents the intra-band Rashba coupling, while 
i^g OI couples bands with different n y indices. In addition to 
SOI, a critical ingredient for realizing Majorana fermions in 
semiconductor nanowires is represented by the Zeeman field. 
We consider that the Zeeman splitting T is generated by ap- 
plying a magnetic field oriented along the wire (i.e., along the 
x-axis), T = gjiBB x /2. The corresponding matrix element 
in the low-energy basis are 



eeman *rr, 



(5) 



where a = —a. 

The third key ingredient is the proximity-induced supercon- 
ductivity (SC). After integrating out the SC degrees of free- 
dom, the SC proximity effect is described by a self-energy 
term that has the form l20l 



(6) 
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FIG. 7. (Color online) Density of states p{u) (top) and end-of-wire 
local density of states p en d (bottom) for a system with single band 
occupancy. The traces (offset for clarity) represent constant magnetic 
field cuts through the color scale diagrams shown in Fig. [6] The DOS 
and LDOS are calculated using Eqns. (HI and (|9|10|>, respectively. 



where 7 is the effective SM-SC coupling, r x and r z are Pauli 
matrices in the Nambu space, A = lmeV is the pair po- 
tential of the bulk SC, and ( is a proximity-induced shift of 
the chemical potential. In the present calculations we take 
£ = 0. As shown in detail in Ref. |20l within the static 
approximation yAo — u2 ~~ ^ me self-energy becomes 
w — 7/A0CJ — JcTyTy and the low-energy physics of 
the SM nanowire with proximity-induced SC can be described 
by an effective Bogoliubov-de Gennes Hamiltonian. This ap- 
proximation is valid, strictly speaking, at energies much lower 
than A , but represents a very good approximation even for 
E ~ Aq/2. Explicitly, the matrix elements of the effective 
BdG Hamiltonian can be written as 



#BdG(™, n') = Z [e n 6 nn ' + Ta x S nn f + (H§ 01 ) nnf ] r z 



A(TyTy, 



(7) 



where e n is given by Eq Q and the SOI matrix elements are 
given by ([?]). Note that energy scale for the SM nanowire is 
renormalized by a factor Z = (I + 7/A0) -1 due to the SC 
proximity effect. This renormalization is determined by the 
term in the self-energy ([6]) that is proportional to u (in the 
static approximation). The pairing term in Eq. ^ is derived 
from the corresponding contribution to the self-energy ([6]) and 
is proportional to the induced pair potential A = 7 A / (7 + 
Ao) = 250/ieV. The effective model described by Eq. ^ is 
solved numerically 

Density of states (DOS) and local density of states (LDOS). 
Diagonalizing the BdG Hamiltonian generates the eigenstates 
®v = (uv,v v ) and the corresponding energies E v with 
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FIG. 8. (Color online) Magnetic field traces corresponding to Fig. 1 
(top and middle panels) from the main text. The dominant features 
above the induced gap A = 250/xeV are due to states from the low- 
energy bands. Note that these features are absent in the single-band 
case (see Fig. [7]). The chemical potential corresponds to the "delo- 
calized" regime, hence in the LDOS the are no signatures associated 
with the closing of the gap at the critical point. 

E_ v = -E v , where v = ±1,±2,... and E 1 < E 2 < 

E — 3< The components u v and v v on the Nambu spinor 

represent the particle and hole contributions, respectively. The 
density of states is defined as 



-Vim 

TT 



E v 



- IT] 



IT] 



(8) 



J2n a u t( n i &)u u (n, a) and a similar expres- 

|2 



where \u v \ 2 

sion holds for \v u \ z . To calculate the LDOS, we define the 
real-space functions u u (i,a) = ^Z n ^n,a(^)u u (n^ a) and , 
similarly, v u (i, a). Using these functions, the LDOS can be 
expressed as 

\u v {i,(j)\ 2 



I E Im 



y<l,, 



E v + irj 



\Vv{i,(T)\ 2 

) + E v + it] 



(9) 

Note that J2 i p(i,uj) = p(u). We note that in this work we are 
not interested in the LDOS at a specific point in the wire de- 
fined on the atomic scale, but rather in the average LDOS near 
the end of the nanowire, which is more closely related to the 
experimentally measured differential tunneling conductance. 
Consequently we define the end-of-wire (average) LDOS as 



(10) 



jZy i x ^iend 



where C is a normalization constant and z e nd = ^end/&, 
with a being the lattice constant, defines the characteristic 
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FIG. 9. (Color online) Magnetic field traces in the "localized" 
regime. The bottom panel represents the same quantity as Fig. 4 
from the main text. The peaks that reveal the closing of the gap at 
the TQPT are generated by states from the top band that are localized 
near the end of the wire. 

length scale of the end-of-wire region. In the calculations 
/ en d ~ 20nm, but the results do not change qualitatively if 
this quantity is varied in the range l-10 2 nm. 

The single-band case and the relative magnitude of vari- 
ous features in the DOS and LDOS. Our main findings are: i) 
In a Majorana nanowire with a chemical potential close to the 
bottom of a band (i.e., in the "delocalized" regime) the end-of- 
wire LDOS does not manifest any signature of the gap closing 
though the TQPT. ii) In a multi-band Majorana nanowire, the 
end-of-wire LDOS exhibits strong features at an energy scale 
associated with the induced SC gap A. These features have a 
weak magnetic field dependence and are generated by states 
from the low-energy bands. In the case of single-band occu- 
pancy, n top = 1, property i) still holds, while ii) is not appli- 
cable, as there is no other band with lower energy. To illustrate 
this case, we calculate the DOS and end-of-wire LDOS for a 
system with one occupied band and /i = 0. The results are 
shown in Fig. [6] Notice that the Majorana peak can be clearly 
seen in the LDOS, but there are no features associated with the 
closing of the gap. Also, in contrast with the multi-band case 
(see Fig. 1 in the main text), there are no features associated 
with the induced gap A. 

If we are interested in the relative strength of various fea- 
tures in the DOS and LDOS, it is rather difficult to obtain a 
good estimate based solely on the color scale diagrams. In- 
stead, it is more convenient to calculate p(uu) and p e nd(^) 
trace at different magnetic fields. For completeness, we show 
these traces for three relevant parameter regimes: a) the sin- 
gle occupancy, "delocalized" regime is shown in Fig. [7] [to 
be compared with color scale plot shown in Fig. |6), b) the 
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FIG. 10. (Color online) Spatial profiles of the three lowest energy 
states in a disordered nanowire. The blue (dark gray) lines corre- 
spond to the lowest energy state. Note that in the topologically triv- 
ial phase (top panel) the lowest energy state is localized in some re- 
gion inside the wire, while for a Zeeman field corresponding to the 
topological SC phase (middle) there are two clear peaks at the ends 
of the wire corresponding to the Majorana bound states. Increasing 
the strength of the disorder potential (bottom) results in effectively 
breaking the wire into smaller segments with Majorana bound states 
localized at their ends. 

multi-band, "delocalized" regime is shown in Fig. [8] [to be 
compared with Fig. 1 from the main text], and c) the multi- 
band, "localized" regime is shown in Fig. [9] [to be compared 
with Fig. 4 from the main text]. 

We point out that the LDOS is closely related to the experi- 
mentally measurable tunneling conductance, but not necessar- 
ily proportional to it in all situations. However, it is difficult 
to think of cases where the LDOS manifests no contribution 
from low energy states (i.e., it is characterized by gap), while 
the tunneling conductance reflects the presence of certain in- 
gap features. As such, our main qualitative finding in this 
work, i.e. that the Majorana mode may show up in the ex- 
periment without the obvious closure of the SC gap, remains 
valid independent of whether the theory focuses on the LDOS 
or the tunneling conductance. 

Role of disorder. Next, we show that the presence of disor- 
der in a system with parameters corresponding to the delocal- 
ized regime does not modify our conclusion. In other words, 
if the gap closing signature is strongly suppressed in the end- 
of-wire LDOS of a clean wire, adding disorder will not gen- 
erate such a signature. The reason for this behavior stems 
from the properties of the low-energy states in the presence 
of a disorder potential [20]. For a typical disorder realiza- 
tion, the low-energy states that are extended in a clean system 
become localized within random segments of the nanowire. 
Nonetheless, localization near the end of the wire is accidental 
and, typically, the contribution of these states to the end-of- 



FIG. 11. (Color online) Dependence of the density of states (top 
panel) and the end-of-wire LDOS (bottom) on the Zeeman field for 
a disordered wire with the same parameters as in Fig. [10| (L x = 
3/xm, A/i = 0, and Vdis = Vdisi). Similar to the clean case, the 
closing of the quasiparticle gap is clearly seen the the total DOS, but 
it is strongly suppressed in the LDOS. The weakly Zeeman field- 
dependent features that are dominant in the LDOS are due to states 
from the lower-energy bands that have a large amplitude at the end 
of the wire. 

wire LDOS remains negligible compared to the contribution 
of the Majorana bound state or the contributions from states 
associated with low-energy bands, which have large ampli- 
tudes near the end of the wire in both clean and disordered 
wires. To illustrate this behavior, we consider a wire of length 
L x = 3/im and a chemical potential at the bottom of the forth 
band (A/i = 0) in the presence of a random potential Vdi s (*) 
with a characteristic length scale [ 20] of the order of the wire 
width, L y . The dependence of the amplitude of the three low- 
est energy states on the position along the wire is shown in 
Fig. 
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Note that the amplitude of the low-energy states go 
smoothly to zero at the ends of the wire, except the Majorana 
bounds state (blue line in the middle panel) that is character- 
ized by a sharp rise. 

The spatial properties of the low-energy states translate di- 
rectly into the dependence of the end-of-wire LDOS on the 



Zeeman field (see Fig. 11 and Fig. 12). Similar to the clean 
case, the strongest contributions to the end-of-wire LDOS 
come from states from the lower energy bands that disperse 
weakly with the Zeeman field. If a low-energy state from the 
Majorana band (the top occupied band) is accidentally local- 
ized near the end of the wire, it will will generate a contribu- 
tion to the end-of-wire LDOS that disperses with the mag- 
netic field. Such contributions can be identified in Figures 1 1 
and [12] However, in order to have a clear gap closing sig- 
nature (e.g., similar to that shown in the lower panel of Fig. 
[9]), it is necessary that the lowest energy mode be localized 
near the end of the wire in the topologically trivial case. This 
requirement is not satisfied for a typical disorder realization. 
Hence, for a system in the "delocalized" regime, the end-of- 
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FIG. 12. (Color online) Zeeman field traces for the DOS (top) and 
the LDOS (bottom) shown in Fig. [TT] The traces are offset for clar- 
ity. Due to the presence of disorder, a few states with energy lower 
than the induced gap Ai n d = 250/ieV are localized near the end of 
the wire and generate features in the LDOS. However, similar to the 
clean case, there is no gap-closing signature. 



wire LDOS will be characterized by the absence of a gap clos- 
ing signature, even in the presence of disorder. 



Role of finite temperature. Finite temperature does not af- 
fect our conclusion adversely at all, in fact, reinforcing our 
finding. More specifically, finite temperature produces broad- 
ening effects, but does not shift spectral weight. Conse- 
quently, any feature in the measured tunneling conductance 
that could be associated with the closing of the gap is sharpest 
at T = 0, while finite temperature reduces or suppresses it. 
Nonetheless, if the weight of a gap-closing feature is compa- 
rable to that of the Majorana ZBCP, it cannot be completely 
suppressed by temperature broadening without destroying the 
ZBCP itself. Furthermore, a feature that is not present at 
T = is unlikely to appear at finite temperature. Thus, finite 
temperature cannot provide a mechanism for the non-closure 
of the gap at the topological quantum phase transition, which 
is the mystery being addressed in our work. It should be em- 
phasized that the actual temperature regime used in the exper- 
iments is substantially below the SC gap value, hence it has a 
very small quantitative impact on the results presented in this 
work. 



